A Bayesian inversion for emissions and export productivity across the end-Cretaceous boundary.

The end-Cretaceous mass extinction was marked by both the Chicxulub impact and the ongoing emplacement of the Deccan Traps flood basalt province. To understand the mechanism of extinction, we must disentangle the timing, duration, and intensity of volcanic and meteoritic environmental forcings. In this study, we used a parallel Markov chain Monte Carlo approach to invert for carbon dioxide (CO2) and sulfur dioxide (SO2) emissions, export productivity, and remineralization from 67 to 65 million years ago using the LOSCAR (Long-term Ocean-atmosphere-Sediment CArbon cycle Reservoir) model. Our results closely match observed and proxy data and suggest decoupled CO2 and SO2 emissions, a two-step decline in export productivity with a protracted recovery, and no clear volatile impulse at the boundary. More broadly, our methods provide a potential path forward for efficient parallel inversion of complex Earth system models.

Continental flood basalts drive Phanerozoic extinctions.

Refinements of the geological timescale driven by the increasing precision and accuracy of radiometric dating have revealed an apparent correlation between large igneous provinces (LIPs) and intervals of Phanerozoic faunal turnover that has been much discussed at a qualitative level. However, the extent to which such correlations are likely to occur by chance has yet to be quantitatively tested, and other kill mechanisms have been suggested for many mass extinctions. Here, we show that the degree of temporal correlation between continental LIPs and faunal turnover in the Phanerozoic is unlikely to occur by chance, suggesting a causal relationship linking extinctions and continental flood basalts. The relationship is stronger for LIPs with higher estimated eruptive rates and for stage boundaries with higher extinction magnitudes. This suggests LIP magma degassing as a primary kill mechanism for mass extinctions and other intervals of faunal turnover, which may be related to [Formula: see text], Cl, and F release. Our results suggest continental LIPs as a major, direct driver of extinctions throughout the Phanerozoic.

Thermochronologic constraints on the origin of the Great Unconformity.

The origin of the phenomenon known as the Great Unconformity has been a fundamental yet unresolved problem in the geosciences for over a century. Recent hypotheses advocate either global continental exhumation averaging 3 to 5 km during Cryogenian (717 to 635 Ma) snowball Earth glaciations or, alternatively, diachronous episodic exhumation throughout the Neoproterozoic (1,000 to 540 Ma) due to plate tectonic reorganization from supercontinent assembly and breakup. To test these hypotheses, the temporal patterns of Neoproterozoic thermal histories were evaluated for four North American locations using previously published medium- to low-temperature thermochronology and geologic information. We present inverse time-temperature simulations within a Bayesian modeling framework that record a consistent signal of relatively rapid, high-magnitude cooling of ∼120 to 200 °C interpreted as erosional exhumation of upper crustal basement during the Cryogenian. These models imply widespread, synchronous cooling consistent with at least ∼3 to 5 km of unroofing during snowball Earth glaciations, but also demonstrate that plate tectonic drivers, with the potential to cause both exhumation and burial, may have significantly influenced the thermal history in regions that were undergoing deformation concomitant with glaciation. In the cratonic interior, however, glaciation remains the only plausible mechanism that satisfies the required timing, magnitude, and broad spatial pattern of continental erosion revealed by our thermochronological inversions. To obtain a full picture of the extent and synchroneity of such erosional exhumation, studies on stable cratonic crust below the Great Unconformity must be repeated on all continents.

Heterogeneous Hadean crust with ambient mantle affinity recorded in detrital zircons of the Green Sandstone Bed, South Africa.

The nature of Earth's earliest crust and the processes by which it formed remain major issues in Precambrian geology. Due to the absence of a rock record older than ∼4.02 Ga, the only direct record of the Hadean is from rare detrital zircon and that largely from a single area: the Jack Hills and Mount Narryer region of Western Australia. Here, we report on the geochemistry of Hadean detrital zircons as old as 4.15 Ga from the newly discovered Green Sandstone Bed in the Barberton greenstone belt, South Africa. We demonstrate that the U-Nb-Sc-Yb systematics of the majority of these Hadean zircons show a mantle affinity as seen in zircon from modern plume-type mantle environments and do not resemble zircon from modern continental or oceanic arcs. The zircon trace element compositions furthermore suggest magma compositions ranging from higher temperature, primitive to lower temperature, and more evolved tonalite-trondhjemite-granodiorite (TTG)-like magmas that experienced some reworking of hydrated crust. We propose that the Hadean parental magmas of the Green Sandstone Bed zircons formed from remelting of mafic, mantle-derived crust that experienced some hydrous input during melting but not from the processes seen in modern arc magmatism.

Constraining crustal silica on ancient Earth.

Accurately quantifying the composition of continental crust on Hadean and Archean Earth is critical to our understanding of the physiography, tectonics, and climate of our planet at the dawn of life. One longstanding paradigm involves the growth of a relatively mafic planetary crust over the first 1 to 2 billion years of Earth history, implying a lack of modern plate tectonics and a paucity of subaerial crust, and consequently lacking an efficient mechanism to regulate climate. Others have proposed a more uniformitarian view in which Archean and Hadean continents were only slightly more mafic than at present. Apart from complications in assessing early crustal composition introduced by crustal preservation and sampling biases, effects such as the secular cooling of Earth's mantle and the biologically driven oxidation of Earth's atmosphere have not been fully investigated. We find that the former complicates efforts to infer crustal silica from compatible or incompatible element abundances, while the latter undermines estimates of crustal silica content inferred from terrigenous sediments. Accounting for these complications, we find that the data are most parsimoniously explained by a model with nearly constant crustal silica since at least the early Archean.

U-Pb constraints on pulsed eruption of the Deccan Traps across the end-Cretaceous mass extinction.

Temporal correlation between some continental flood basalt eruptions and mass extinctions has been proposed to indicate causality, with eruptive volatile release driving environmental degradation and extinction. We tested this model for the Deccan Traps flood basalt province, which, along with the Chicxulub bolide impact, is implicated in the Cretaceous-Paleogene (K-Pg) extinction approximately 66 million years ago. We estimated Deccan eruption rates with uranium-lead (U-Pb) zircon geochronology and resolved four high-volume eruptive periods. According to this model, maximum eruption rates occurred before and after the K-Pg extinction, with one such pulse initiating tens of thousands of years prior to both the bolide impact and extinction. These findings support extinction models that incorporate both catastrophic events as drivers of environmental deterioration associated with the K-Pg extinction and its aftermath.

Neoproterozoic glacial origin of the Great Unconformity.

The Great Unconformity, a profound gap in Earth's stratigraphic record often evident below the base of the Cambrian system, has remained among the most enigmatic field observations in Earth science for over a century. While long associated directly or indirectly with the occurrence of the earliest complex animal fossils, a conclusive explanation for the formation and global extent of the Great Unconformity has remained elusive. Here we show that the Great Unconformity is associated with a set of large global oxygen and hafnium isotope excursions in magmatic zircon that suggest a late Neoproterozoic crustal erosion and sediment subduction event of unprecedented scale. These excursions, the Great Unconformity, preservational irregularities in the terrestrial bolide impact record, and the first-order pattern of Phanerozoic sedimentation can together be explained by spatially heterogeneous Neoproterozoic glacial erosion totaling a global average of 3-5 vertical kilometers, along with the subsequent thermal and isostatic consequences of this erosion for global continental freeboard.

Potassic, high-silica Hadean crust.

Understanding Hadean (>4 Ga) Earth requires knowledge of its crust. The composition of the crust and volatiles migrating through it directly influence the makeup of the atmosphere, the composition of seawater, and nutrient availability. Despite its importance, there is little known and less agreed upon regarding the nature of the Hadean crust. By analyzing the 87Sr/86Sr ratio of apatite inclusions in Archean zircons from Nuvvuagittuq, Canada, we show that its protolith had formed a high (>1) Rb/Sr ratio reservoir by at least 4.2 Ga. This result implies that the early crust had a broad range of igneous rocks, extending from mafic to highly silicic compositions.

A record of deep-ocean dissolved O2 from the oxidation state of iron in submarine basalts.

The oxygenation of the deep ocean in the geological past has been associated with a rise in the partial pressure of atmospheric molecular oxygen (O2) to near-present levels and the emergence of modern marine biogeochemical cycles. It has also been linked to the origination and diversification of early animals. It is generally thought that the deep ocean was largely anoxic from about 2,500 to 800 million years ago, with estimates of the occurrence of deep-ocean oxygenation and the linked increase in the partial pressure of atmospheric oxygen to levels sufficient for this oxygenation ranging from about 800 to 400 million years ago. Deep-ocean dissolved oxygen concentrations over this interval are typically estimated using geochemical signatures preserved in ancient continental shelf or slope sediments, which only indirectly reflect the geochemical state of the deep ocean. Here we present a record that more directly reflects deep-ocean oxygen concentrations, based on the ratio of Fe3+ to total Fe in hydrothermally altered basalts formed in ocean basins. Our data allow for quantitative estimates of deep-ocean dissolved oxygen concentrations from 3.5 billion years ago to 14 million years ago and suggest that deep-ocean oxygenation occurred in the Phanerozoic (541 million years ago to the present) and potentially not until the late Palaeozoic (less than 420 million years ago).

Volcanic-plutonic parity and the differentiation of the continental crust.

The continental crust is central to the biological and geological history of Earth. However, crustal heterogeneity has prevented a thorough geochemical comparison of its primary igneous building blocks-volcanic and plutonic rocks-and the processes by which they differentiate to felsic compositions. Our analysis of a comprehensive global data set of volcanic and plutonic whole-rock geochemistry shows that differentiation trends from primitive basaltic to felsic compositions for volcanic versus plutonic samples are generally indistinguishable in subduction-zone settings, but are divergent in continental rifts. Offsets in major- and trace-element differentiation patterns in rift settings suggest higher water content in plutonic magmas and reduced eruptibility of hydrous silicate magmas relative to dry rift volcanics. In both tectonic settings, our results indicate that fractional crystallization, rather than crustal melting, is predominantly responsible for the production of intermediate and felsic magmas, emphasizing the role of mafic cumulates as a residue of crustal differentiation.

Statistical geochemistry reveals disruption in secular lithospheric evolution about 2.5 Gyr ago.

The Earth has cooled over the past 4.5 billion years (Gyr) as a result of surface heat loss and declining radiogenic heat production. Igneous geochemistry has been used to understand how changing heat flux influenced Archaean geodynamics, but records of systematic geochemical evolution are complicated by heterogeneity of the rock record and uncertainties regarding selection and preservation bias. Here we apply statistical sampling techniques to a geochemical database of about 70,000 samples from the continental igneous rock record to produce a comprehensive record of secular geochemical evolution throughout Earth history. Consistent with secular mantle cooling, compatible and incompatible elements in basalts record gradually decreasing mantle melt fraction through time. Superimposed on this gradual evolution is a pervasive geochemical discontinuity occurring about 2.5 Gyr ago, involving substantial decreases in mantle melt fraction in basalts, and in indicators of deep crustal melting and fractionation, such as Na/K, Eu/Eu* (europium anomaly) and La/Yb ratios in felsic rocks. Along with an increase in preserved crustal thickness across the Archaean/Proterozoic boundary, these data are consistent with a model in which high-degree Archaean mantle melting produced a thick, mafic lower crust and consequent deep crustal delamination and melting--leading to abundant tonalite-trondhjemite-granodiorite magmatism and a thin preserved Archaean crust. The coincidence of the observed changes in geochemistry and crustal thickness with stepwise atmospheric oxidation at the end of the Archaean eon provides a significant temporal link between deep Earth geochemical processes and the rise of atmospheric oxygen on the Earth.